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Several recent works on quantum criticality beyond the Landau- Ginzburg- Wilson paradigm have 
led to a number of field theories, potentially important for certain two-dimensional magnetic in- 
sulating systems, where criticality is not very well understood. This situation highlights the need 
for non-perturbative information about criticality in two spatial dimensions (three space-time di- 
mensions), which is a longstanding challenge. As a step toward addressing these issues, we present 
evidence that the 0(4) vector model is dual to a theory of Dirac fermions coupled to both SU(2) 
and U(l) gauge fields. Both field theories arise as low-energy, long-wavelength descriptions of a 
frustrated XY model on the triangular lattice. Abelian boson-vortex duality of the lattice model, 
together with the emergence of larger non-abelian symmetry at low energies, leads to this rare ex- 
ample of duality in two spatial dimensions involving non-abelian global symmetry and fermions, 
but without supersymmetry. The duality can also be viewed as a bosonization of the Dirac fermion 
gauge theory. 
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I. INTRODUCTION 

Recent theoretical progress on quantum criticality in 
two spatial dimensions (d = 2) has elucidated a num- 
ber of critical phenomena beyond the Landau-Ginzburg- 
Wilson (LGW) paradigm, particularly in the context of 
d = 2 magnetic insulators i 1 ' 2 ' 3 ' 4 ' In the LGW paradigm, 
one identifies one or more order parameters and uses 
them to construct an effective field theory based on sym- 
metry alone, and studies criticality in this theory. The 
field theories that arise outside the LGW framework 
make up a much broader class than those arising within 
it, often involve fermions and/or gauge fields, and in gen- 
eral are less understood than those theories arising within 
the LGW framework. 

Most analytical understanding of quantum criticality 
in d = 2 is based on perturbative calculations in large-n 
and d = 3 — e expansions, and it would be extremely use- 
ful to obtain non-perturbative information that reduces 
our reliance on these formal limits. This is even true 
within the LGW framework, although, for example, the 
0(n) model is relatively well understood after extensive 
numerical and analytical studies. 5 The field theories that 
arise beyond the LGW paradigm have been much less 
studied, and non-perturbative information is more ur- 
gently needed. Toward this end, in this paper I argue that 
the 0(4) model (with some anisotropy terms) is identical 
at low energies to an apparently quite different theory 
involving fermions and gauge fields, which is closely re- 
lated to field theories describing one of the main classes 
of non-LGW criticality, namely critical spin liquids. This 
relation is a rare example of duality in D = d + 1 = 3 
space-time dimensions involving non-abelian global sym- 
metries, non-abelian gauge fields, and fermions, without 
supersymmetry. 

The proposed relation, which is the main result of this 
paper, is that the 0(4) vector model is dual to a the- 
ory of Nf = 2 four-component Dirac fermions coupled 
to both SU(2) and U(l) gauge fields. We dub the lat- 



ter theory QCED3, as it is in a sense a combination of 
D — 3 quantum electrodynamics (QED3) and quantum 
chromodynamics (QCD3) with N c = 2 colors. To obtain 
this relation, we begin with an XY antiferromagnet on 
the triangular lattice. This model can be treated directly 
by a Landau theory approach, and this leads to the 0(4) 
field theory, with various anisotropy terms breaking the 
0(4) symmetry down to that of the original model. Al- 
ternatively, using abelian boson- vortex duality)^ the XY 
model is mapped onto a dual lattice model, from which 
we obtain QCED3 as a low-energy effective description. 
The non-abelian duality between the 0(4) model and 
QCED3 thus emerges at low energy as a descendant of 
abelian boson-vortex duality on the lattice. 

This result builds in a crucial way on work of Alicea 
et. al., where it was proposed, beginning with the same 
frustrated XY model, that the 0(4) vector model is dual 
to Nf — 2 two-component Dirac fermions coupled to 
a U(l) Maxwell gauge field and a U(l) Chern-Simons 
gauge fields Later, this result was obtained from a differ- 
ent point of view by Senthil and Fisher, who also argued 
that the 0(4) model with a topological term (0-term at 
6 — 7r) is dual to the same U(l) gauge theory but with no 
Chern-Simons term^ In these examples, either the origi- 
nal or dual theory involves a topological term, which pre- 
vents the development of a controlled large-n or d = 3 — e 
expansion. Because such expansions are possible for the 
dual partners without topological terms, these results are 
very useful for understanding their partners that do have 
topological terms. However, because only one dual part- 
ner can be directly analyzed, such dualities do not lead 
to much additional understanding of the critical behav- 
ior itself. The result presented here differs in the crucial 
respect that both the original and dual theories lack topo- 
logical terms, and thus both admit controlled expansions, 
potentially allowing greater insight into the critical prop- 
erties. 

In general, two field theories are dual when: (1) the 
two theories have identical low-energy physics, and (2) 
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the variables in the two theories are mutually nonlocal. 
That is, the fields representing one theory cannot be writ- 
ten as a local function of the fields in the other, and vice 
versa. The latter condition excludes trivial, local changes 
of variables. The classic example of such a duality for 
D = 3 is between the XY model and the abelian Higgs 
model. In general, it is known how to construct explicit 
duality relations for models with abelian global symme- 
tries and abelian gauge fields - without fermions - in ar- 
bitrary space-time dimension^ However, as in the field 
theories of interest here, many of the most interesting 
cases involve non-abelian global symmetries, non-abelian 
gauge fields and/or fermions, where no general techniques 
exist to construct useful duality relations. The primary 
exceptions are in D = 2, where bosonization 1 "^ can be 
thought of as a duality, and in a number of supersym- 
metric field theories in higher dimensions . 12 ' 13 ! 14 ' 15 The 
results of this paper are a step toward understanding du- 
ality in a less restrictive context. 

It is important at this stage to make a distinction be- 
tween duality of field theories and duality of critical fixed 
points. By duality of field theories, we shall mean that as 
the parameters of the two theories are varied, they have 
the same phases and critical points, and all the same 
low-energy, long-wavelength behavior. Duality of fixed 
points, on the other hand, is a stronger statement, and is 
most easily defined by discussing its meaning in the con- 
text of the 0(4)-QCED3 duality. The 0(4) model has 
a Wilson-Fisher fixed point which can be accessed in a 
controlled fashion by a large-n expansion, where one con- 
siders the 0(n) model. Also, QCED3 has a critical fixed 
point that can be accessed in the large-iVy expansion; 
we assume that this critical point survives down to the 
value Nf = 2 of interest here. By duality of fixed points, 
we would mean in this case that the 0(4) Wilson-Fisher 
fixed point is identical to the QCED3 fixed point. How- 
ever, it must be emphasized that this statement need not 
hold for the two field theories to be dual. It could be the 
case that the 0(4) and QCED3 fixed points are distinct, 
and, in each of the two field theories, both fixed points 
are present. The 0(4)-QCED3 duality as proposed is to 
be understood primarily as a duality between field theo- 
ries. Whether it also leads to a duality between the 0(4) 
and QCED3 fixed points is a more challenging question, 
and the available evidence does not allow us to reach a 
firm conclusion. However, it is plausible that the two 
fixed points are identical. 

We now briefly outline the remainder of the paper. In 
Sec. Hi] we introduce the frustrated triangular lattice XY 
model, and review the treatment of Ref. |8| up to the point 
where it deviates from our approach. In particular, we 
discuss a direct Landau theory approach, which leads to 
the 0(4) field theory (Sec. Ill Aj) . and the dual lattice 
model resulting from the boson-vortex duality transfor- 
mation (Sec. Ill B|) . Next, we describe our route from the 
dual model to a continuum field theory (Sec. IIII A|) . and 
proceed to give a more precise statement of the proposed 
duality (Sec. IIII B[) . In the following two sections, we 



present the principal evidence for the duality: In Sec. IIVI 
we set up a dictionary between operators in the 0(4) and 
QCED3 field theories, and in Sec.|V]we demonstrate that 
the three stable phases of the original XY model can be 
realized in QCED3. We conclude in Sec. I VII with a dis- 
cussion of open issues and possible directions for future 
work. Three appendices contain various technical details. 



II. PRELIMINARIES AND PRECURSORS 

A. Direct analysis of XY model 

The initial steps needed to obtain the 0(4)-QCED3 
duality have already been carried out by Alicea et. al., 8 
so we review their results here and in Sec. Ill Bl while in- 
troducing notation that will be important later on. At 
a certain point, which we identify in Sec. Ill Bl our treat- 
ment deviates from theirs. As we shall very briefly out- 
line, following their route leads instead to the duality 
between the 0(4) model and Dirac fermions coupled to 
both Maxwell and Chern-Simons U(l) gauge fields. Our 
route to the 0(4)-QCED3 duality is described beginning 
in Sec. iHIAl 

The starting point, here and in Ref. |g, is an XY anti- 
ferromagnet on the triangular lattice, whose Hamiltonian 
is 

H = Uj2 n l + J J2 ™*(4>r ~ 4>r>)- (1) 

r {rr') 

Here, J > and the second sum is over nearest-neighbor 
bonds. On every site of the triangular lattice r there is 
a U(l) quantum rotor, with angular position e*^ r and 
integer- valued angular momentum n r , which satisfy the 
commutation relation [n, e 1 ^] = e 1 ^ . This can be thought 
of as an effective theory for a S — 1 XY antiferromag- 
net. Alternatively, by the usual quantum-classical corre- 
spondence, its partition function describes a classical XY 
antifcrromagnet of stacked triangular lattices, which has 
been extensively studie d 16 ' 17 ' 18 ' 19 (see also Ref. [H for a 
more detailed overview of the literature). In this context, 
it is known that the model can realize three phases (pro- 
vided next- neighbor ferromagnetic exchange is included) , 
which are: (1) paramagnet (2) coplanar 120° magnetic 
order (Fig. [T]) and (3) two different collinear magnetic or- 
ders (Fig.[T]and Ref.|8|). The paramagnet, coplanar state, 
and one of the two collinear states meet at a multicriti- 
cal point (upon tuning both J and the further-neighbor 
exchange), which is the Wilson-Fisher fixed point of the 
0(4) vector model. Either collinear state can arise ad- 
jacent to the multicritical point, depending on the sign 
of a dangerously irrelevant 6th-order term in the Landau 
theory (discussed below) , and the distinction between the 
two collinear states will not play an important role in the 
present discussion. 

The symmetries of Eq. ([1} will play a crucial role in 
our analysis, so we enumerate them here. The triangular 
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FIG. 1: Depiction of the coplanar 120° ordered state (a), and 
one of the two collinear states (b) . The filled circles represent 
sites with zero ordered moment. The other collinear state is 
depicted in Ref. @. 



lattice space group is generated by a translation by one 
lattice constant in the a;-direction (T x ), a counterclock- 
wise rotation by 7r/3 about a lattice site and a 
reflection lZ y : (x.y) — > (x, — y). The operators n r and 
<f) r transform as scalars under the space group. There is 
also a charge-conjugation (or spin flip) symmetry 



C : (j) r — » — 4> r . 
And the antiunitary time-reversal operation 



T 
T 



(2) 
(3) 



(4) 
(5) 



Finally there is the U(l) phase rotation, which sends 
(j) r — > (f> r + a. To distinguish it from other symmetry 
groups that will arise later on, we shall refer to this sym- 
metry as U(1)xy- 

To expose the 0(4) structure, one goes to a path 
integral description of the partition function in terms 
of (j> r (r), where r is imaginary time, and focuses on 
those configurations of the phase field with lowest action. 



These are given by e 



icf) r 



where Q = (47r/3)a; (the lattice constant is set to unity), 
and zi t 2 are slowly varying functions of r and r. The 
vectors ±Q lie at the corners of the hexagonal Brillouin 
zone and are the ordering wavevectors of the 120° state. 
The action on #x,2 of the microscopic symmetries is enu- 
merated in Appendix [B] 

Constrained by the microscopic symmetries, one then 
writes a continuum effective Lagrange density 



C = \d„Z\ 2 + rZ^Z + u(Z^Z) 2 + v 4 \ Zl \ 2 \ 



(6) 



where terms of order |z| 6 and higher have been discarded, 
and Z T = (z\Z2). For v 4 = 0, the model's continuous 
symmetry is SO(4), which is broken down to U(l) x U(l) 
for V4 7^ 0. The 6th order terms break the continuous 
symmetry down to \J(1)xy- The v 4 and r terms are rele- 
vant perturbations to the 0(4) multicritical point, while 
the allowed 6th order (and higher) terms are irrelevant. 
At the mean-field level, the coplanar state arises when 
r < and v 4 > 0, and has (zx) 7^ and (z 2 ) = (or vice- 
versa). The collinear state obtains for V4 < and has 



I (zi) I = I (z 2 ) I 7^ 0. Depending on the relative phase of z\ 
and z 2 two different collinear ordering patterns are possi- 
ble; this phase is determined by the sign of a dangerously 
irrelevant 6th order term. 

Some of the properties of the 0(4) critical point are 
known from numerical simulations and high-order per- 
turbative calculations; see Ref. HI for an extensive dis- 
cussion (and tabulation) of critical properties of O(n) 
models. Approximately, the critical exponent v m 0.75, 
and 77 w 0.027. These exponents can be translated into 
scaling dimensions of fields by dimZ = (1 + r))/2 ps 0.51 
and dirnZ^Z = 3- 1/u w 1.67. 

It will be convenient to use the fact that SO (4) ~ 
[SU(2) x SU(2)]/Z 2 . We thus introduce the matrix 



z x -z 2 

Z2 Zx 



(7) 



A general SO(4) rotation is realized by Z — * UlZUr, 
where Ul.r <= SU(2). We shall thus refer to left and right 
SU(2) rotations, denoted by SU(2)x i ij, with conserved 
currents J^' R - In terms of the fields, 



jR 



Z^(d^Z) 



(8) 
(9) 



where cr = (a 1 , a 2 , cr 3 ) is a vector of the 2x2 Pauli matri- 
ces. The XJ(1)xy symmetry is a subgroup of SU(2)a and 
is generated by 2{Jq ) z , while Jq generates rotations of 
the two slowly varying fields Zx and z^ into one another. 

One of the major pieces of evidence for the 0(4)- 
QCED3 duality will be a dictionary identifying operators 
in the two field theories. To that end, we now enumerate 
some important operators in the 0(4) model. Operators 
can be labeled by (£l>£r), the total angular momentum 
quantum numbers of SU(2)j, x SU(2)#. It will some- 
times be useful to also specify (rriL,rnR), the projection 
quantum numbers of the angular momenta along the z- 
axis. Note that 2m r is the U(1)xy charge. We have 
already discussed the boson field Z itself, which trans- 
forms as (1/2, 1/2), and the currents and Jjf, which 
transform as (1,0) and (0, 1), respectively, and are both 
vectors under Lorentz rotations (i.e. space-time rota- 
tions). We shall also consider the following 9 Lorentz 
scalar bilinears in Z: 



N 
I 
I 



Z^crZ 

Z T {ia 2 )crZ 

(!)*■ 



(10) 

(11) 
(12) 

= 1 multiplet. We 

= Z^Z and O v = 

\zx\ 2 \z 2 \ 2 - (l/6)(ZtZ) 2 . O r has £ L = l R = 0, and O v is 
a member of a £l — £r — 2 multiplet with tol = mji = 0. 
These play an important role as they are the terms in the 
Lagrangian that must be tuned to reach the 0(4) critical 
point. Next, defining the 0(4) vector <^,- (i = 1, . . . , 4) by 



Together these make up a £l — £r 
shall also consider the operators O r 
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FIG. 2: Relation between the dual honeycomb lattice (dashed 
lines) and the triangular lattice of the original XY model 
(solid lines) . Honeycomb lattice sites correspond to triangular 
lattice plaquettes. 



21 = 4>i + i<$>2 and z-i = 4>3 + i(f>4, we define the operator 
£ = itijki£ t j.v\4>td ll 4>jd v 4>kd\(j)i, which is odd under T and 
lZ y . fn a nonlinear sigma model version of the 0(4) field 
theory, where the constraint Z'Z = 1 is imposed, if £ 
is integrated over space-time and added to the action it 
becomes the topological 6*-term. 20 While € does not have 
any special topological significance as a local operator, 
we shall refer to it as the topological density. 



B. Dual model 

In Ref. la, boson-vortex duality was applied to the 
Hamiltonian Eq. (J). This is a straightforward proce- 
dure for such a lattice XY model; here we review some 
crucial aspects, and more details are found in Ref.|H. The 
degrees of freedom of the dual model are vortices, which 
are bosons residing on the sites of the dual honeycomb 
lattice (Fig. [5]), and the vector potential and electric field 
of a non-compact U(l) gauge field, which reside on the 
nearest-neighbor bonds of the honeycomb lattice. Cru- 
cially, the vortices are at half-filling (an average of one- 
half vortex per honeycomb lattice site). This is a direct 
consequence of the frustration in the original XY model: 
the XY exchange term of Eq. |T]) can be rewritten as 
— J^2(rr') cos^r^c^r'+Arr'), where the flux associated 
with the non-fluctuating vector potential A rr ' is half a 
flux quantum, or tt, for each triangular plaquette. This 
background flux has the consequence of forcing in half a 
vortex per site (on average) in the dual theory. 

For technical convenience, we shall take the vortices 
to experience a hardcore repulsion, and only allow zero 
or one vortex on each honeycomb site. Such a modifica- 
tion of microscopic parameters will not affect the univer- 
sal low-energy properties of phases and critical points in 



which we shall be interested. 
The dual Hamiltonian is 

Hduai - (/'^[(Vxa)o] 2 + J'^e^ 

O (ff') 

- t v J2 [e w -X+ V + H.c] . (13) 

{ff') 

Here, f labels the sites of the dual honeycomb lattice, e^/ 
and aff> are, respectively, the electric field and vector po- 
tential residing on the honeycomb link joining f and f'. 
On the same link these satisfy the canonical commuta- 
tion relation [e, a] — i. The first term is a sum over all 
hexagonal plaquettes, and (V x a)o is the lattice curl, 
which is the discrete line integral of a ff i taken around 
the perimeter of a given hexagonal plaquette. The lat- 
ter two terms are sums over nearest-neighbor honeycomb 
links. Because the vortices are hardcore bosons, their cre- 
ation, destruction and number operators are represented 
using the S — 1/2 spin operators (vf,vf,v^), satisfying 
the usual commutation relation [u£,t> f ,] = iSpp^^Vp. A 
vortex at site f is created by vf = vf + iv% and de- 
stroyed by Vf = (Vf )'. The vortex number operator is 
N v (f) = vf+1/2. The Hamiltonian is supplemented with 
the Gauss' law constraint 

(dive) P = v z f , (14) 

where (div e)? — X)f'~r e "" ^ s * ne lattice divergence of 
e f ft. The parameters of the dual Hamiltonian are U' ~ 
U, J' ~ J, and the vortex hopping t v . It should be 
emphasized that Eq. (|13|) is not an exact rewriting of the 
original XY model, and instead should be thought of as 
a low-energy effective theory. Conversely, the XY model 
of Eq. ([I]) can also be thought of as a low-energy effective 
theory for the dual model Eq. (fT3"|) . 

It is important to spell out the connection between 
the U(l) gauge field and the original XY model degrees 
of freedom. The magnetic flux represents the V(1)xy 
charge density, 

-!-(V x a) ~ n r , (15) 

Z7T 

where O is the honeycomb hexagon surrounding the tri- 
angular lattice site r. Also, the electric field is related to 
the 11(1) xy current by 

sin(27re F r') ~ sin(</> r — 4> r '), (16) 

where (r,r') is the unique triangular lattice bond cross- 
ing the honeycomb link (f, f'). In the dual representa- 
tion, XJ(1)xy charge conservation is represented as the 
conservation of magnetic flux. Moreover, insertion of 
U(1)xy charge (as by acting with e 1 ^) is represented 
by the insertion of a quantized 27r magnetic flux. Such 
flux insertions are space-time magnetic monopole events, 
and will play a crucial role in our analysis. 

Under the space group symmetry, e^f and aff trans- 
form as pseudovectors, and as a pseudoscalar. vf 
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transforms as a scalar under T x and and lZ y : v f — > 

, where r' is the image of f under the reflection. Under 
charge-conjugation both e and a are odd, while 



C 
C : 



-v? 



(17) 
(18) 



Finally, under time-reversal a is odd while e is even, and 



T 
T : 



(19) 
(20) 



The challenge at this stage is to use the dual Hamilto- 
nian to construct a continuum low-energy effective the- 
ory, which would then provide a dual description of the 
0(4) critical point. A natural approach would be to begin 
with Eq. (|13[) and construct a functional integral in terms 
of the phase field of the vortices v + ~ e lB . However, as 
the vortices are at half-filling, the resulting functional in- 
tegral is plagued with Berry phase terms, and does not 
directly lead to a useful continuum limit. Instead, it is 
necessary to choose some other set of variables that is 
amenable to a continuum description. 

Alicea et. al. chose to represent the vortices in terms of 
fcrmions using statistical transmutation, with the price of 
introducing a Chern-Simons gauge field ctrf'-— In a flux- 
smearing mean-field treatment, the vortices are fermions 
at half-filling on the honeycomb lattice (with zero back- 
ground magnetic flux), and thus have a massless Dirac 
dispersion. To include fluctuations about the mean-field 
state, one re-couples the Dirac fermions to both the 
Chern-Simons gauge field and the Maxwell gauge field 
a that arose in the boson-vortex duality. The resulting 
Euclidean Lagrangian is 



-T-OLpfLavxdvOlX 
47T 



(21) 



The first term is the kinetic energy of the fermions and 
contains the coupling to the two gauge fields. The second 
term dictates that the gauge field obeys Maxwell dy- 
namics, and the last term is the Chern-Simons term for 
a M . The coefficient of the Chern-Simons term is precisely 
that needed to attach 2ir flux and transmute fermions 
into bosons (and vice versa) . The ellipsis represents other 
(important) perturbations consistent with the underlying 
microscopic symmetries. 

While this field theory is an intriguing dual representa- 
tion of the 0(4) model, the presence of the Chern-Simons 
term seriously hinders any direct analysis of it. It was 
conjectured in Ref. H that the Chern-Simons term can 
simply be dropped without affecting the critical proper- 
ties, but the arguments in favor of this conjecture are not 
conclusive, and are questionable given the later results of 
Ref. Motivated in part by a desire to avoid these is- 
sues, here we pursue a different approach that also leads 
to a fermionic gauge theory representation of the 0(4) 
model, but without any topological terms. 



III. DUALITY BETWEEN 0(4) MODEL AND 
QCED3. 



A. Vortex Fractionalization Route to Dual 
Effective Theory 

The present approach also begins with the dual Hamil- 
tonian Eq. (flU)) . The challenge is to construct a contin- 
uum effective theory that deals with the half-filling of the 
vortices, but also avoids the difficulties associated with 
statistical transmutation and the Chern-Simons term. To 
do this, we look to the theory of spin liquids in S = 1/2 
Heisenberg models, which are after all also models of half- 
filled hardcore bosons. One route to describe spin liquids, 
in particular critical spin liquids that lack a spin gap, is to 
formally represent hardcore bosons (or spins) as bilinears 
of fermionic slave particles^ Therefore it is reasonable 
to hope that a similar approach can describe criticality 
- in particular, the 0(4) critical point - in the present 
model, and we shall argue that this is indeed the case4£ 

We represent the vortices using the fermionic operators 
fr a , where a — 1, 2: 



— ff 1 ff2 



(/rl/fl ffzffz)- 



(22) 
(23) 



With the local constraint fl a f fa = 1 this change of vari- 
ables provides an exact rewriting of the model. In these 
variables there is a local SU(2) gauge redundanc y 22 ^ 3 
which can be exposed by defining 



= f fk 



(24) 



where Cf = 1 (—1) for r in the A (B) sublattice. The 
fermions satisfy the local constraint equations iplf^ipf = 
0, where (i — 1,2,3) are the 2x2 Pauli matrices. 
This is simply the condition that the SU(2) gauge charge 
is zero at every lattice site. The vortex operators and 
Vf can be written in manifestly SU(2) gauge-invariant 
forms, and so all physical operators are gauge invariant. 

To proceed, we pass to a functional integral represen- 



tation, where the action is S = S- 



U(l) 



Sf, where 



5*u(i) = JJjJ dT ^2 [ d Tdff' ~ (Aa )rf'] 2 

{fr') 

+ U' J drJ2(^ x a )o- ( 25 ) 



Here, we have introduced the notation (A/)^' = f{r') — 
f(f) for a lattice derivative. The fermionic part of the 
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action is 

S f = / <lr 



4 [dr + 2«o(r, r) + °^ yP 



+ U 



/*-£[■ 



e ia "'u±^ +H.c. 



{rr') 



Vv 



(26) 



Note that we have made the gauge transformation Off — > 
Off' + 7r to change the sign of the last term. We have also 
introduced the Lagrange-multiplier held a l (i = 1,2,3), 
which enforces the SU(2) gauge constraint and can be 
thought of as the time-component of the SU(2) gauge 
field. 

The last term of Eq. (|26p is quartic in the fermion op- 
erators, and to arrive at a candidate low-energy effective 
theory it can be decoupled using a Hubbard-Stratonovich 
field residing on the bonds of the lattice. One searches for 
mean-field saddle points of this field, and each such sad- 
dle point (upon including often important fluctuations) 
leads to a low-energy effective theory. Many distinct ef- 
fective theories can be generated in this way^i and one of 
the challenges of this approach is to decide which theory 
(if any) accurately captures the physics of the model at 
hand. In the present case we are guided by the require- 
ment that the effective theory should be able to repro- 
duce the phases and critical points that are known to be 
present from the analysis of the original XY model. 

Rather than carry out the above mean-field procedure 
explicitly, we shall simply guess the form of the low- 
energy theory. To do this we write down an effective 
lattice gauge theory that reduces to the above model in 
a particular limit. This is equivalent to choosing a par- 
ticular mean-field saddle point and then including the 
fluctuations about it. The effective lattice theory is ob- 
tained by replacing Sf with 



dr 



2 / dT ^2 a o( r ' T ) 

f 

tJdrJ2 [e~ ia ^' /2 4fUffi> f - + H -' 



(27) 



{rr') 



Here Ufr> is the spatial part of the SU(2) gauge field. 
It should be noted that, because the gauge field dyy/ IS 
noncompact, it is perfectly legitimate to have the ob- 
ject e mfF '/ 2 appearing in the action. We also include a 
Maxwell action S g for the SU(2) gauge field, with over- 
all strength proportional to the coupling constant l/.g 2 . 
When g — ► oo and t is small, Uff 1 can be integrated 
out perturbativcly in t. At leading order (i 2 ), one re- 
covers the original dual Hamiltonian given by the action 
S = S*u(i) + Sf, with an additional nearest-neighbor re- 
pulsive interaction between vortices. This interaction is 



not important for our purposes, as it does not change the 
symmetry of the model and is not expected to affect its 
universality class. 

On the other hand, if we take g small and U' large, then 
the fluctuations of both gauge fields are suppressed, and 
in an appropriate gauge Uff ~ 1 and a l rs a??' ~ &o ~ 0. 
In this mean-field limit the fermions are described by the 
Hamiltonian 



Hmft = t ^2 [4 i'f + H - c -] i 



(28) 



(rr') 



which is simply nearest-neighbor hopping of half-filled 
fermions on the honeycomb lattice. The corresponding 
low-energy theory is given by focusing on the excitations 
near the Dirac nodes and reintroducing the coupling to 
the gauge fields. The resulting Lagrangian density is 



£qced3 



•J 



$ (29) 



Here we have introduced the continuum Dirac field $, 
which is related to Vv as discussed in Appendix [A] $ 
is an eight-component object; these eight components 
arise from the two-component nature of ijif, the two 
bands needed to represent each Dirac node, and the two- 
component flavor index corresponding to the two distinct 
nodes in the Brillouin zone. It is useful to define three 
different sets of Pauli matrices acting in this 8-component 
space; each set corresponds to its own type of SU(2) ro- 
tations. The /j, 1 Pauli matrices act in the SU(2) gauge 
space and generate gauge transformations. The t 1 Pauli 
matrices act in the band index, or Lorentz, space, and 
generate Lorentz transformations. Finally, the v l Pauli 
matrices act in the SU(2) flavor space and generate fla- 
vor rotations. $ resides in the tensor product of these 
three SU(2) spaces, and products of different types of 
Pauli matrices (which commute) should be understood 
as matrix tensor products. The action of the various 
Pauli matrices on <& is given explicitly in Appendix [A] It 
is convenient to think of $ as composed of Nf = 2 fla- 
vors of four-component Dirac fermions, where each flavor 
transforms as a doublet under SU(2) gauge rotation and 
Lorentz transformations. Both flavors carry the same 
U(l) charge of 1/2 under the dual gauge field a M . 

We have also introduced the field strength of the SU(2) 
gauge field, 



(30) 



and the matrices 7^ are defined in terms of the t % Pauli 
matrices as 



7 ^(r 3 ,r 2 ,-r 1 ). 
Finally, we have defined 



(31) 



(32) 
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A number of theories similar to £qced3, also involving 
Dirac fermions coupled to gauge fields, are interesting in 
the context of non-LGW criticality. This interest stems 
from the fact that the microscopic symmetries can (in 
some cases) be enough to forbid the addition of any rel- 
evant perturbations to the action, and the fixed point 
thus describes a stable critical phase ^£21 Such criti- 
cal phases have been discussed in a variety of physical 
settings 25.26,27,28,29.30,31,32,33,34 

The global symmetries of £qced3 play a crucial role 
in our discussion. Aside from D = 3 Poincare invariance, 
the continuous global symmetry is SU(2) x U(l). The 
global SU(2) consists of rotations between the two flavors 
of Dirac fermions generated by the v % Pauli matrices. We 
thus dub it SU(2)ir, and it has the conserved current 

Jl = §7 M i/$. (33) 

The global U(l) is simply the U(1)xy symmetry. Its 
realization in £qced3 follows directly from the boson- 
vortex duality transformation, and in particular the iden- 
tification of a M magnetic flux and U(l)xy charge density 
[Eq. (p~5|) ] . The associated conserved current is thus 

ja = 7r£nv\d v a\, (34) 

the flux of the a M gauge field. 

The field theory £qced3, like other D = 3 theories of 
massless Dirac fermions coupled to gauge fields, is solv- 
able in the large- Y/ limit ^LS^SML When Nf -> 00 
the fluctuations of the gauge fields are suppressed, and 
for most purposes the physics is identical to that of non- 
interacting fermions. Expanding about this limit, cor- 
relation functions can be calculated order-by-order in 
1/Nf, and it is found that various operators acquire 
anomalous dimensions. The principal result, then, is that 
the large-iV/ expansion describes an interacting critical 
fixed point which we assume survives down to the 
case of interest, Nf = 2. Below, we shall refer to the 
Nf = 2 incarnation of the large-AT/ fixed point as the 
QCED3 fixed point, which should be distinguished, of 
course, from the QCED3 field theory. 

B. Statement of the proposed duality 

A brief statement of the proposed duality is that the 
0(4) model and QCED3 possess identical low-energy 
physics. Below, we shall elaborate on the meaning of 
this statement in order to give a more precise statement 
of the duality. We then outline the approach underlying 
the evidence for the duality, which is described in the 
following sections. 

It is useful to remark that both the 0(4) model 
and QCED3, provided arbitrary perturbations consistent 
with the underlying microscopic symmetries are added 
to each, are expected to be valid low-energy effective de- 
scriptions of the original model. This is expected simply 



because both were derived from the same microscopic 
starting point. This means, in particular, that any phase 
or critical point of cither the direct or dual low-energy 
theory is a phase or critical point that presumably ex- 
ists somewhere in the parameter space of the original 
XY model. It is not immediately obvious, however, that 
the 0(4) model and QCED3 describe phases and critical 
points in the same part of parameter space. This is the 
crucial fact that needs to be established for the duality 
to hold. 

With this remark in mind, the proposed duality can 
be precisely stated as follows: Beginning with £qced3, 
by adding operators consistent with the underlying mi- 
croscopic symmetries one can tune the theory to an 0(4) 
critical point identical to the Wilson-Fisher fixed point 
of the 0(4) model. We shall call the resulting fine-tuned 
Lagrangian £qced3- Moreover, the effective Lagrangian 

£ cff = C c QCEm + rO r + v6 v , (35) 

where O r and O v are QCED3 operators identified be- 
low in Sec. IIV1 is identical at low energy to the 0(4) 
model Lagrangian of Eq. ©. These statements essen- 
tially amount to asserting that, indeed, the 0(4) model 
and QCED3 describe the same region of the original 
model's parameter space. 

Using the terminology introduced in Sec. U this is a 
duality between field theories. It is natural to ask if there 
is also a duality between the 0(4) Wilson-Fisher fixed 
point and the QCED3 critical point. That is, are the 
0(4) and QCED3 fixed points identical or distinct? We 
do not have a definitive answer to this question, but we 
will see it is plausible that the two fixed points are dual. 
This is discussed further in Sec. lVIl Most of the evidence 
presented in Secs. HVI andfVlpertains to the duality of field 
theories. However, we do present some information on 
scaling dimensions from the large- Nf expansion for the 
QCED3 fixed point. These results are not relevant for 
establishing the duality between field theories, but they 
do provide some information about a potential duality of 
fixed points. 

In Sec. HVI we set up a dictionary of operators between 
QCED3 and the 0(4) model. In order to do this, it 
will be helpful to exploit the large- Nf understanding of 
the QCED3 fixed point; this provides a controlled un- 
derstanding of the field content of QCED3, which is our 
primary concern in establishing the duality between field 
theories. Next, in Sec. [V] we describe how to access the 
stable phases of the original XY model in the QCED3 
description. 



IV. EVIDENCE FOR THE DUALITY: 
OPERATOR DICTIONARY 

To construct a dictionary between QCED3 and 0(4) 
model operators, we begin by identifying the continuous 
global symmetries of the two field theories. We identify 
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the SU(2) F flavor symmetry of QCED3 with SU(2) L of 
the 0(4) model, which means that the currents and 
should also be identified. Next, we identify jj* with 
2{J^) Z , as these are the realizations of the U(1)xy con- 
served current in QCED3 and the 0(4) model, respec- 
tively. 

As QCED3 does not have manifest SU(2) x SU(2) 
symmetry, it is less clear how to identify the remain- 
ing components of Jj}. An important part of the pro- 
posed duality is that, at the 0(4) critical point, the 
SU(2) L x U(1)xy global symmetry of QCED3 is enlarged 
to SU(2)i x SU(2) fl . So, in particular, U(1)xy is en- 
larged to SU(2)ij, which happens quite explicitly in the 
0(4) model. Moreover, it should not be surprising that 
it is possible to enlarge SU(2) x U(l) to SU(2) x SU(2) 
by a suitable tuning of parameters, since both groups at 
least have the same Cartan subalgebra. We shall identify 
candidate QCED3 partners of the remaining components 
of jj} below. 

To continue setting up our dictionary, it will be con- 
venient to break QCED3 operators into two classes. The 
first class consists of all operators carrying zero U(1)xy 
charge, or, equivalently, zero magnetic flux of the dual 
gauge field et M . We refer to operators in the first class 
as non-monopole operators. The second class contains 
all monopole operators, which do carry an flux. The 
flux is quantized in multiples of 2n, simply because the 
U(1)xy charge is quantized in the original XY model. 
Non-monopole operators can be easily represented in 
terms of the fermion and gauge fields. On the other hand, 
as is typical for topological disorder operators, monopole 
operators are more difficult to represent in terms of the 
fields of the theory. However, they can be constructed us- 
ing the state-operator correspondence of conformal field 
theory^ and we shall take advantage of this approach 
here. 

To be identified, two operators certainly must trans- 
form identically under Lorentz transformations, SU(2)x, 
(or, equivalently, SU(2)f), and U(1)xy- For non- 
monopole operators, in each case we have also verified, 
using the results enumerated in Appendix [B] that each 
pair of identified operators transforms identically under 
all the microscopic symmetries. For monopole operators, 
on the other hand, it is only known how to partially de- 
termine the action of microscopic symmetries, up to a few 
unknown parameters i 8 ' 39 i 40 Here we shall do this, follow- 
ing the approach of Ref. 40|. In each case the transforma- 
tions of the corresponding 0(4) model operator can be 
determined completely using the results of Appendix [Bl 
and are consistent with the partial results for the QCED3 
counterpart. 

We now proceed to identify some of the important 0(4) 
model operators with QCED3 counterparts. 

0(4) bilinear AT = Z^crZ. This is perhaps the sim- 
plest operator to identify with a QCED3 partner. We 
identify N with 

N = -i$v<P. (36) 



Both N and TV transform as a vector under SU(2)^, a 
scalar under Lorentz rotation, and carry zero U(1)xy 
charge. 

In the limit Nf — > oo, the scaling dimension A^ of AT 
approaches 2. An inspection of the diagrams involved in 
the 1/Nf correction shows that the contributions from 
the U(l) and SU(2) gauge fields-'- 7 simply add together, 
and 

1 28 

A * = 2 -3^ + ° (1/A ^ (3?) 

While this result is not to be trusted quantitatively, 
the qualitative trend that < 2 is believed to be 

reliable^i^i This should be compared to the fact that, 
in the large-n expansion of the 0(n) model, Ajv = 
I + 0(l/n). It is reasonable that these scaling dimen- 
sions become the same for the case of interest (n = 4 and 
Nf = 2). 

0(4) field Z. Since Z carries U(1)xy charge of unity, 
it is a monopole operator carrying 2tt flux. To construct 
27r-monopoles in QCED3, we need to make use of the 
state-operator correspondence following Ref. [H. Acces- 
sible treatments of the state-operator correspondence can 
be found in Refs.l42l andl43l. Briefly, for a D — 3 Lorentz 
and scale invariant theory, such as the QCED3 fixed 
point, the state-operator correspondence states that local 
operators of the field theory are in one-to-one correspon- 
dence with quantum states of the same theory quantized 
on the two-dimensional unit sphere. Scaling operators 
(i.e. eigenoperators of a renormalization group transfor- 
mation) are mapped to eigenstates of the Hamiltonian on 
the 2-sphere, and the scaling dimension of the operator 
is equal to the energy of the corresponding state. 

Here we are primarily concerned with using the state- 
operator correspondence to construct monopole opera- 
tors, by constructing the corresponding eigenstates of 
the Hamiltonian on the 2-sphcre. This can be done for 
QCED3 in the large- AT/ limit, where the fluctuations of 
both gauge fields are suppressed and the Nf — ► oo Hamil- 
tonian is simply that of Nt flavors of non-interacting 4- 
component fermions on the unit sphere. As Nt is reduced 
from infinity down to the case of interest Nf — 2, gauge 
fluctuations will modify the states and change their ener- 
gies (and hence their scaling dimensions). However, even 
for Nf — 2, the states constructed with gauge fluctua- 
tions suppressed still correspond to local operators (al- 
though no longer to scaling operators), and we still use 
them to understand the field content of Nf = 2 monopole 
operators. 

We shall see below that it will be useful to think in 
terms of 7V/2 = 2Nf two-component fermions, each of 
which is a Lorentz doublet. A monopole operator with 
flux 2wq corresponds to a state of these fermions with 
a background flux of 2nq on the sphere. Since the 
fermions carry et M charge of 1/2, they feel only a flux of 
Ttq. 

Let us now consider a 27r-flux monopole operator, 
where the fermions feel a flux of 7r from a M ; this violates 
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the Dirac quantization condition and is thus apparently 
problematic. To illustrate the problem, recall that the 
monopole's gauge field can be represented using a Dirac 
string carrying the 27r-flux away from the center of the 
monopole in an infinitesimally thin solenoid. However, 
since the fermions feel only a 7r-flux from the solenoid, it 
is a physical object, and the resulting object is thus not 
a local operator. One might try to eliminate the Dirac 
string using the mathematical technology of sections and 
the monopole harmonics of Wu and Yang^ but Dirac's 
quantization condition still enters there as a requirement 
that the transition function be single valued, and the 
problem is not avoided. 

However, the Dirac quantization condition can be re- 
paired in the present case if flux in the SU(2) gauge 
field is also present, and so this must happen in order 
to get a local operator with unit U(1)xy charge. In 
particular, we can consider putting a 27r-flux monopole 
in - half of the 2-component fermions feel this as 
7r flux, and the other half feel it as —tt flux. (Other 
strengths of monopoles can also be considered, but these 
either lead to states that will have higher energy - and 
thus higher scaling dimension - or continue to violate 
the Dirac quantization condition.) Therefore, combined 
with the overall 7r-flux coming from a M , there will be 
Nf 2-component fermions moving in a background 27r 
flux, and Nf fermions moving in zero flux. An issue that 
immediately arises is that the resulting state is not in- 
variant even under global SU(2) gauge transformations. 
We return to this below, after discussing the structure of 
the state in the fixed gauge where there is a monopole 
in a?, and a}, = a 2 ., = 0. It should be noted that this 
configuration of the SU(2) gauge field is not itself topo- 
logically stable, but is induced by the topologically stable 
monopole in the U(l) gauge field. 

To understand the structure of the 27r-monopole states 
and thus the corresponding operators, we shall need the 
spectrum of a single 2-component Dirac fermion on the 
2-sphere, moving in a uniform background flux 27r/. The 
energy levels for |/| > 1 are^ 

E p = ±v/p 2 +p|/|, (38) 

where p is a nonnegative integer. For / = the form of 
the spectrum is the same, but p is restricted to be positive 
and thus there are no E p — states. These levels are 
(2j p + l)-fold degenerate, where j p = (1/2)(|/| — 1) + p, 
and transform in the (2j p + l)-dimensional representation 
of the SU(2) Lorentz group (i.e. rotation symmetry of 
the 2-sphere). 

We now specialize to Nf = 2, the case of interest. We 
represent $ in terms of iV/2 = 4 two-component fermions 
T] aa , where a, a = 1,2 and $ T = (77^, 77^, 77^, ?yj 2 ). The 
a index transforms as a doublet of SU(2)l, and the a 
index as a doublet of SU(2) gauge rotations. This means 
that the a — 1 (a = 2) fermions carry charge of 1/2 
(—1/2). Then 7711 and 7721 feel flux / = 1, while 7712 and 
7722 feel / = 0. The spectrum is as illustrated in Fig. [31 
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FIG. 3: Illustration of the 27r-monopole state whose corre- 
sponding operator is the QCED3 partner of the 0(4) field 
Z. The spectrum of each 2-component fermion ri aa is shown. 
The fermions 7711 and 7721 both feel a net 27r flux through the 
sphere. The energy level spectrum for each of these fermions 
is symmetric about zero energy, where for each fermion there 
is a single zero energy state. The 7712 and 7722 fermions feel no 
flux, and also have a spectrum symmetric about zero energy, 
but with no zero energy state. All negative energy levels are 
filled, and one of the two zero mode states is filled to obtain a 
state with zero total charge under both the a M and gauge 
fields. 



In order for the state to correspond to a gauge- 
invariant operator, it must carry zero a M charge. More- 
over, we shall see below that it must also carry zero 
charge. Now, both a M and a^, and hence the correspond- 
ing charges, change sign under the charge conjugation 
symmetry C, which acts as a particle-hole transforma- 
tion on the 77-fermions. Crucially, C acts trivially on the 
a index of r\ aa (see Appendix [Bj) , so it does not exchange 
fermions feeling flux with those feeling no flux. Acting 
with C on the state illustrated in Fig. [3] results in another 
state of the same schematic form - the only change is in 
the occupation of the zero mode states - which must 
therefore have the same charges as the original state. 
The only consistent possibility is that both a M and 
charges are zero, as needed. Other states with the cor- 
rect charges, and higher energy, can be obtained from 
this state by moving fermions between levels. 

Now we return to the question of gauge invariance of 
the 27r-flux monopole state. To correspond to a gauge- 
invariant local operator, the state should be invariant 
under general U(l) and SU(2) gauge transformations. 
We can construct a gauge-invariant state starting with 
the gauge-fixed state discussed above, denoted by l^o) ■ 
(Which of the two zero modes is filled is not important 
for this discussion, so we just focus on a single state.) 
The gauge-invariant state is constructed by integrating 
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over all possible gauge transformations: 



[dGi] / [dG 2 ]G x G 2 \i> ), 



(39) 



where G\ and G2 are unitary operators implementing 
U(l) and SU(2) gauge transformations, respectively, and 
the integrals are taken over all such unitary transforma- 
tions. This state is clearly gauge invariant, but we 
need to check that it does not vanish. This is done in Ap- 
pendix[Cl where it is shown that zero a M and charge is 
a necessary and sufficient condition for \tp) not to vanish. 

We are now in a position to discuss the quantum num- 
bers of the 27r-fiux monopole states depicted in Fig. El 
and thus those of the corresponding operator. The neg- 
ative energy single-particle levels, when fully occupied, 
are invariant under both Lorentz rotation and SU(2)i, 
flavor. However, the fermions feeling / = 1 have 2 
zero-energy Lorentz-singlet levels, which are filled with 
a single fermion. This implies that the overall state 
transforms as a Lorentz singlet, but as a doublet under 
SU(2)l. We thus denote the corresponding monopole 
operator as the two-component object Z, where Z T — 
(zi,Z2), and Z transforms as a doublet under SU(2)l 
and a singlet under Lorentz rotation. The complex con- 
jugate Z* carries U(1)xy charge —1 and corresponds to 
a state on the sphere with — 2ir flux in a M . These are 
the same quantum numbers carried by the 0(4) field Z 
(and its conjugate Z*), and it is therefore reasonable to 
identify Z and Z. 

To work out the transformations of Z under micro- 
scopic symmetries, we follow the procedure of Ref. 
(which is based on that of Ref. d) , using the transforma- 
tions for the fermion field <!> enumerated in Appendix |5] 
The basic idea is to use the facts that Z is a SU(2) l dou- 
blet and carries unit V(1)xy charge to write down the 
most general transformation laws possible, which are 





Z - 


■» e^ T exp 


ir/3 


Z - 


-> e l ' t,R (j 1 Z 


Ky 


Z - 


-> e i<ht -Z 


c 


Z - 


-> e^^Z* 


T 


z - 


c r Z\ 



(40) 

/?, .-. : Z — , '•'< : n A Z (41) 

(42) 
(43) 
(44) 

where <Pt, (/>r, 4>tz and <j) c are arbitrary phases, and cq- — 
±1. We can partially determine these parameters by de- 
manding that the action of the space group on Z sat- 
isfy algebraic relations satisfied by its generators, and by 
making redefinitions of Z. The relations TZ 2 = R®/ 3 = 1 



>R 



= 0, 7r and 



Also, T X R 2 ,„T X R 

1 x 7T/J x 



R ~ 
-2 
t/3 



tttir/3, where 
Z' = e^e^Z, 



give the constraints 
n R = 0, 1, ... ,5 

gives 4>t — 0. Next, we redefine Z 
where a = —<f) c /2 ) which has the effect of setting <fi c = 
in the above transformation laws. If hr < 3, we choose 
(3 = 0, while for Hr > 3 we choose j3 = 7r/2, which has 
the effect of sending riR — > tir — 3 in the transformation 



laws above. To summarize, the most general transforma- 
tion is then characterized by 4>t — § c — 0; 4>ti = 0: n i 
ct = ±1; and 4>r = Tmn/3 where tir = 0,1,2. If we 
choose cr = — 1 and 4>n = tir = 0, we obtain the same 
transformation laws as for the 0(4) field Z, so it is in- 
deed consistent with microscopic symmetries to identify 
Z and Z . 

Bilinears J = Z T (ia 2 cr)Z and I*. These objects 
have V(1)xy charge of ±2. The QCED3 partners can 
thus be represented as states on the sphere where 
has 47r flux. In such a background gauge field, each r\ aa 
fermion feels In flux, and no flux of the SU(2) gauge 
field is needed to satisfy the Dirac quantization condi- 
tion. There are now four zero-energy single particle lev- 
els, and the lowest energy monopole state is obtained 
by filling all negative energy levels and filling the zero- 
energy levels with two fermions. We denote the vacuum 
state for ±47r flux with no zero-energy levels filled by |±), 
and let c\ a create a fermion in the zero-energy level cor- 
responding to i] aa . Focusing on +4ir flux, we consider 
the class of states c^ a cL|+). There are a total of 6 such 
states. Three of these form a triplet under global SU(2) 
gauge transformations, and thus do not correspond to a 
local operator. The other 3 states are a singlet under 
global SU(2) gauge transformations, and a triplet under 
SU(2)i. These states can be labeled by \P) (i = 1,2,3) 
and written 



(icr 2 ) a0 [a l (ia 2 )} ab cl a cl 



+ ■ 



(45) 



We denote the operator corresponding to these states by 
the vector J, while its complex conjugate I* is repre- 
sented by the corresponding states for flux of —An. 

Just as for I (I*), I (J*) is a vector under SU(2) L , a 
Lorentz singlet, and carries U(1)xy charge of +2 (—2). 
Transformations of I under microscopic symmetries can 
be worked out following the same procedure as above for 
Z, and can be chosen to agree with those of the 0(4) 
model operator I. Alternatively, if we think of J as 
the composite operator / ~ Z T (i<r 2 cr)Z, then J and J 
transform identically, simply following from the identical 
transformations of Z and Z. 

Conserved current Jj}. We have already identified 
2{jj}) z with the gauge flux , but it remains to identify 
the other components [Jj}) x , [Jj}) y ■ We define linear 
combinations 



(46) 



so that (J a) carries U(l)xr charge ±2 and corresponds 
to a QCED3 monopole operator. 

To construct the corresponding monopole operator us- 
ing the state-operator correspondence, we let the op- 
erators d\ a ^ p create a fermion in an energy E p level 
for p = ±1. Here the index fi transforms in the 
triplet representation of the Lorentz group. The states 
e p.v\d aa ^ + a' l ab d baX _ \P) have the correct quantum num- 
bers to be identified with {Jj}) + ■ However, there is an- 
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other natural set of candidate states. Denning the gauge- 
triplet states 

|G 4 ) = {io 2 ) ab [a\io 2 )Uci a cl p \+), (47) 

we see that the states e t _ LV xd\ au+ a l af3 d al}x _\G l ) also have 
the correct quantum numbers to be identified with 
(J^) + . While it is not clear which of these operators 
should be identified with (J^) + , the more important 
point is that we have found at least one operator with 
the correct quantum numbers. 

0(4) "mass" term O r = Z^Z. This object must be 
identified with a QCED3 operator that is invariant under 
all symmetries of the field theory. Moreover, at the 0(4) 
fixed point, it is the most relevant such operator. Using 
the QCED3 fixed point as a guide, the most relevant such 
operators in the large-iV/- limit (all with scaling dimen- 
sion 4) are the Maxwell terms for the two gauge fields, 
and the following set of three quartic terms: ($$) 2 , 
(<|>i/<i>) 2 and (i>/x$) 2 . These quartic terms are indepen- 
dent, but to form a complete basis for all singlet quartic 
terms the additional operator (3>7 M $) 2 must also be in- 
cluded. However, it can be shown that this operator ac- 
tually has scaling dimension 6 in the Nf — > oo limit.— As 
Nf is reduced from infinity, some of the above dimension- 
4 operators are expected to become more relevant. In- 
deed, to identify the QCED3 and 0(4) fixed points, one 
of them must lower its dimension to about 1.67, the di- 
mension of Z'Z at the 0(4) critical point, while the oth- 
ers must remain irrelevant (dimensions greater than 3). 
To have a clearer indication whether this is the case, it 
would be useful to calculate l/Nf corrections to the di- 
mensions of these operators. While similar calculations 
have appeared in the literature before^SiiS Nf = oo shifts 
in scaling dimension [analogous to the fact that ($7 M $) 2 
has dimension 6 here] are present but were not taken into 
account; this will modify the results. We hope that fu- 
ture work will resolve this issue, which would be useful in 
a variety of situations where similar field theories arise. 

Anisotropy term O v = \z 1 \ 2 \z 2 \ 2 - (1/6)&Z. This 
operator is a Lorentz singlet, carries zero U(1)xy charge, 
and belongs to a 4 = 2 multiplet of SU(2)i (with 
n%L = 0), so, as with O r it is natural to identify it with 
a QCED3 operator quartic in the fermion fields. We 
have classified quartic terms according to their transfor- 
mations under the Lorentz group and SU(2)i, and find 
that there are two independent terms with the correct 
quantum numbers: 

Ro = (#^ 3 $) 2 - \($v<P) 2 (48) 

Ri = (^ 3 7m*) 2 - i(^7^) 2 - (49) 

The 0(4) operator O v should be identified with the more 
relevant of Rq and R\ . 

Topological density €. We identify the 0(4) oper- 
ator £ with the QCED3 bilinear £ = — which is 



also odd under spatial reflections and time-reversal, but 
invariant under other symmetries. 

It may seem that this identification spells trouble for a 
potential duality between 0(4) and QCED3 fixed points. 
£ is quartic in the 0(4) field and has three derivatives, 
thus appearing strongly irrelevant, while £, as a fermion 
bilinear, naively appears likely to be relevant. However, 
two points are in order. First, the 0(4) model operator £ 
cannot be generalized in a natural way to the 0(n) case 
(or to d = 3), so the above intuition about its scaling 
dimension is suspect. Second, the scaling dimension of £ 
to leading order in the large- A^ expansion is 

25fi 

dim£ = 2+ 3^V7 +0(1/ ^ 2) ' (50) 

which, upon inspection of the diagrams involved, can be 
obtained by simply adding the corresponding anomalous 
dimensions for QED3^> and QCD3.— Therefore, gauge 
fluctuations make £ substantially more irrelevant - the 
coefficient of l/Nf is quite large. The striking contrast 
from the behavior of the N bilinears, which become more 
relevant with decreasing Nf, is interpreted physically in 
Ref. [4l|. It is then plausible that £ continues to be irrel- 
evant, perhaps strongly so, down to Nf = 2. 

V. EVIDENCE FOR THE DUALITY: STABLE 
PHASES 

We shall now show that the three stable phases of the 
0(4) model Eq. (J6j) - the paramagnet, and coplanar and 
collinear ordered states - are realized simply in QCED3. 
While we do not have enough control over QCED3 at 
Nf — 2 to determine the phase diagram, the fact that all 
the known phases of the 0(4) model can be realized is an 
important piece of evidence in support of the proposed 
duality. 

We begin by considering the paramagnetic phase, 
which can be represented in the dual description by con- 
densation of vortices^ We need to find a vortex conden- 
sate that does not break any microscopic symmetries, 
and to that end we consider 

U = $ T (^ 2 )(w 2 )(ir 2 )$, (51) 

which creates a single vortex. V is invariant under T x 
and T and is a scalar under R n /^. Under !Z y and C, 
V —> V'. Working in a fixed gauge, condensation of V 
means that (V) = e l6 \V\ ^ 0. If 9 = 0, then clearly 
the condensate preserves the microscopic symmetries. In 
fact, even if 9 ^ 0, the apparent breaking of 1Z V and C 
is a gauge artifact, as these symmetries can be restored 
by supplementing them with an appropriate gauge trans- 
formation. Alternatively, one can always choose a gauge 
such that 9 — 0. 

We can understand the effect of this vortex condensa- 
tion on the fermion spectrum, at the mean-field level, by 
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adding the term (A/2) J d 2 r(V + V^) to the mean-field 
Hamiltonian. It is easily seen that the dispersion rela- 
tion becomes E(k) = ±^/k 2 + A 2 , and a gap of 2A is 
opened. In this gapped state, the SU(2) gauge field will 
become confining, thus removing single fermion excita- 
tions from the spectrum. Also, the U(l) gauge field will 
be in a Higgs phase, and its photon mode and excita- 
tions carrying nonzero XJ(1)xy charge will acquire a gap. 
Therefore all excitations are gapped, no symmetries are 
broken, and we have a description of the paramagnet. 

To access the magnetically ordered phases, we consider 
the Lagrangian C = £qced3 + ot ■ N. At the mean-field 
level (and also at the QCED3 fixed point), the addition 
of a ■ N will open a gap in the fermion spectrum. As a 
result the SU(2) gauge field will become confining, and 
monopole operators of the U(l) gauge field will acquire an 
expectation value - the latter phenomenon corresponds 
to the spontaneous breaking of V(1)xy symmetry and 
XY magnetic order. The photon of the U(l) gauge field 
remains gapless, and corresponds to the spin wave mode 
of the magnetically ordered state. 

To determine the pattern of magnetic order, we note 
that because a ■ N breaks various space group symme- 
tries, it leads to a gapped "vortex insulator," where the 
pattern of lattice symmetry breaking depends on a. The 
vortex insulators that result are the same as those arising 
in the discussion of Ref. [H, so we give only a brief sum- 
mary of the results here. If a = a z z, the resulting state is 
a vortex "charge density wave," where vortices preferen- 
tially occupy one of the two honeycomb sublattices. This 
corresponds to the 120° coplanar magnetically ordered 
state. 8 On the other hand, if a = a x x + a y y, the result- 
ing state is a vortex "valence bond solid" (VBS), where 
vortices hop back and forth preferentially on a subset of 
honeycomb lattice bonds. As detailed in Ref. [H, the vor- 
tex VBS states that arise here correspond to the collincar 
magnetically ordered states. 



VI. DISCUSSION 

We have given evidence for a rare example of non- 
abelian duality (without supersymmetry) in d — 2, be- 
tween the 0(4) model and QCED3 field theories. Both 
were derived as low-energy effective descriptions of the 
same triangular lattice frustrated XY model. The 0(4) 
model arose from a standard Landau theory treatment. 
QCED3 was obtained by combining the abelian boson- 
vortex duality of the lattice model with a fermionic slave- 
particle treatment of the vortex degrees of freedom. By 
setting up an operator dictionary between the two field 
theories, and showing that both can realize the three sta- 
ble phases of the original XY model, we argued that both 
field theories describe the same low-energy sector of the 
XY model and are thus dual. 

As discussed in Sec. IIII B[ this result is to be under- 
stood primarily as a duality between field theories. It 
is natural to ask whether it is also a duality between 



the 0(4) and QCED3 fixed points, assuming that the 
QCED3 fixed point exists for Nf = 2. The two possi- 
bilities are that either the 0(4) and QCED3 fixed points 
are identical, or they are distinct. The latter possibil- 
ity would mean, for example, that the 0(4) model real- 
izes both fixed points somewhere in its parameter space. 
While we are not aware of any evidence for this, it is 
hard to rule out, especially if the Nf = 2 QCED3 fixed 
point is highly unstable and requires tuning of several 
parameters. However, given the available evidence on 
scaling dimensions of operators at the two fixed points, 
it is certainly plausible that they are identical. Numeri- 
cal simulations of QCED3 for various values of Nf could 
potentially shed some light on this issue. 

The 0(4)-QCED3 duality suggests a number of direc- 
tions to pursue a further understanding of d = 2 non- 
abelian duality. One can certainly explore whether sim- 
ilar constructions, perhaps starting from other lattice 
models, lead to other duality relations. Another nat- 
ural question is whether there are connections between 
the present results and d — 2 dualities of supersymmetric 
field theories . 14 ! 15 Perhaps upon suitably breaking super- 
symmetry in a pair of dual theories, the 0(4)-QCED3 
duality - or other, related dualities - can be obtained. 

Finally, it is interesting to remark that the 0(4)- 
QCED3 duality can be viewed as a bosonization of 
QCED3, where the bosonized form is simply the 0(4) 
model. If more examples of similar dualities are discov- 
ered, and better understood, they might eventually be 
useful as a kind of d = 2 bosonization. 
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APPENDIX A: CONTINUUM LIMIT FOR 
FERMIONS 

Consider the mean-field Hamiltonian 

Hmft = t ^ (i>Ur> + H - c -)' ( A1 ) 

(rf') 

where the sum is over nearest-neighbor links of the hon- 
eycomb lattice. We will take the continuum limit for 
the low-energy excitations (which are massless Dirac 
fcrmions), and give the relationship between continuum 
and lattice fields. 

We use the 2-site unit cell labeled by (R,i), with i = 
1,2, and define the Fourier transform 

^ a = -^£ e *-fy Wj (A2) 
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where N c is the number of unit cells in the lattice. The 
Hamiltonian then becomes 



(A3) 



and we have 



d 2 q 



(A15) 



where 

[H(k)] n = [H(k)] 22 

and 

[H(k)] 12 = [H(k)]* 21 = 1 



— ik-a-2 



(A4) 



(A5) 



The Dirac nodes are at k = ±Q, where Q = (4:w/3)x. 
Putting k = ±Q + q and expanding to first order in small 
q we find 



and 



H{Q + q) = -\iQxT 1 +q y r 2 ), 



H(-Q + q) = ^j-(q x T 1 -q y T 2 ), 



(A6) 



(A7) 



where r* are 2x2 Pauli matrices. We define the contin- 
uum spinors 



and 



<?i (?) 



<P2(q) 



^Q+qA 



i>-Q+q.l 
1p-Q+q.2 



(A8) 



(A9) 



Note that (p a (with a = 1, 2) is a 4-component object. 

It is convenient to act on this space with tensor prod- 
ucts of two different kinds of Pauli matrices. The t 1 Pauli 
matrices introduced above act in the space labeling the 
two sites of the unit cell, and the fi z Pauli matrices gen- 
erate SU(2) gauge transformations. For example, the 
action of these Pauli matrices on (f>i{q) is given by 



and 



//^Q+q,2 



la ^Q+q,a 
2a ^Q+q,a 



(A10) 



(All) 



where in each entry there is an implied sum over a = 1,2. 
The continuum Hamiltonian takes the form 



d 2 q 



(2n) 2 



-q x r 



T^ifi+ifiliq^-qyT 2 )^ 



(A12) 

where v is the velocity. To make this look like the con- 
ventional Dirac Hamiltonian we define 



Pi 
92 



(A13) 
(A14) 



We then define the eight-component object 



$ = 



<Pi 
92 



(A16) 



We introduce another set of Pauli matrices, is 1 , acting in 
this two-component flavor space. For example, 



v 1 ® = { ^ 
Vi 



(A17) 



APPENDIX B: SYMMETRIES 

Here we enumerate the action of the microscopic sym- 
metries of the triangular lattice XY model Eq. ([T]) (listed 
in Sec. UTS}, on both the 0(4) model and QCED3 fields. 

The action of the symmetries on the 0(4) model 
field Z T — {z\,z>2) is easily obtained from e l ^ r ~ 
e l Q- r Zl (r,T) + e~ l ®' r z 2 (r, r), the relationship between 
lattice and continuum fields discussed in Sec. Ill Al One 
obtains the following results: 



T x : Z{r, t 



R 



tt/3 

TZy 

c 

T 



Z(r,r 
Z(r,r 
Z(r,r 
Z(r,r 



(- 



exp v 3 

u x Z{r',r) 
Z(r', T ) 
g x Z*(t,t) 
-Z*(r,r). 



Z(r,r) 



(Bl) 

(B2) 
(B3) 
(B4) 
(B5) 



For -Rtt/3 and TZ y , r' is the image of r under the corre- 
sponding symmetry operation. 

To obtain the action of the symmetries on the QCED3 
fermion field <f>, one first obtains their action on the lat- 
tice fermion ipf. The action of each symmetry is cho- 
sen so that (1) the vortex operators u£, which are bilin- 
ears of have the correct transformation properties, 
and (2) Hmft is invariant. Because tpf is not a gauge- 
invariant object, the symmetry group acting on it should 
be thought of as a projective symmetry groupi^i One 
feature of this situation is that, because 7Ymft is invari- 
ant under global SU(2) gauge transformations, with each 
symmetry operation we are free to make an arbitrary 
global SU(2) gauge transformation. The action of the 
symmetries on $ is then easily determined by the rela- 
tions between continuum and lattice fields given in Ap- 
pendix [X] Several analyses of this kind have appeared in 
the literature (see, for example, Ref. [36h . so we simply 
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quote the results: 





/47ri 3 \ 




K/3 ■ *(r) - 




(B7) 


ft, : $(r) - 


■+ -^ 2 )^ 2 )$*(r') 


(B8) 


C : $ - 




(B9) 


T : $ - 


4 (i/l 2 )(iT 2 )(^ 2 )$. 


(BIO) 



Combined with the requirement that the action be in- 
variant, these transformations determine the transforma- 
tions of the SU(2) gauge field a^. We also enumerate the 
transformations of the U(l) gauge field a^. Under T x , 



both and a 



^ are invariant. Under R^/^, both ao and 



a l transform as scalars, while the spatial components 
rotate as vectors. Under the remaining symmetries, 



TZy : (a , ai 


a.2) - 


■* (-oo,-ai,a 2 ) 


(BH) 


TZy : {al,a\, 


4) - 




(B12) 


C 


: a/i - 


■+ -Op 


(B13) 


C ■ (&l,al, 


«£) - 


■* (-ap, ap, -a£) 


(B14) 


T 




■+ -Op 


(B15) 


T 




- a*. 


(B16) 



Note that the transformations of ct\ depend on the arbi- 
trary choice of global SU(2) gauge transformation in the 
corresponding transformation of 



APPENDIX C: NON- VANISHING OF 2tt-FLUX 
MONOPOLE STATE 

Here we show that the gauge- invariant 27r-flux 
monopole state constructed in Eq. (|39| from |-0o) , the 
gauge-fixed state, is nonzero. The Hilbert space of the 
theory on the sphere is a tensor product of fermion and 
gauge field Hilbert spaces, and we can write 



l^> = l/o>®Mfi),oL(n)>, 



(Cl) 



where |/o) gives the state of the fermions, and the lat- 
ter factor the state of the gauge fields in the vector po- 
tential basis, fl = (9, <fi) specifies the angular position 
on the sphere. The vector potential states of the gauge 
field Hilbert space form an orthonormal basis. In spheri- 
cal coordinates, the vector potentials above are given by 
ag = a\ = and 



1 — cos 9 

2sin6» 
1 — cos 9 

2sin6> ' 



(C2) 
(C3) 



This corresponds to a monopole with 2ir flux piercing the 
sphere in both and or„. 

To show that 7^ 0, it is enough to show that 
("00 IV 1 ) 7^ 0. Consider a general gauge transformation 



G = G1G2 acting on \ipo), where Gi is a general U(l) 
gauge transformation and G 2 a general SU(2) gauge 
transformation. Below, we show that either 

G\ih) = (C4) 



G^o) = |/)®|Sp(fi),4,p(^)). (C5) 

where a M a M or a 1 ^ 7^ a* . In the latter case, 
(ipo\G\i/jo) = 0. This implies that 

|V> = J [dG x ] J [dG 2 ]G 1 G 2 |0 o ) = G|0 O ) + W), (C6) 

where G > and (V>o|VO = 0. Therefore {ipoW = C^0, 
and 1^) is indeed nonzero. 

It is clear that \ip) would be zero if either the a M or 
a 3 charge were nonzero. The reason is that the integral 



of Eq. (|C6p is a projection onto gauge-singlet states, if 
either the a M or a£ charge is nonzero, this projection 
must vanish. This shows that it is necessary for these 
charges to vanish in order to have nonzero \ip). 

Let s(f2) be the generator of U(l) gauge transforma- 
tions, and ti(Sl) (i = 1,2,3) the generators of SU(2) 
gauge transformations. Then a general U(l) gauge trans- 
formation can be written G\ = exp(i J*dQA(Q)s(f2)) 
and a general SU(2) gauge transformation is G2 = 
exp(z J dQXi(Cl)ti{fi)). Because |^o) nas zer0 charge un- 
der a M and a 3 , it is invariant under global U(l) gauge 
transformations (A(O) constant), and a global SU(2) 
gauge transformation where Ai = A 2 = and A3 (11) is 
constant. 

It is easy to see that any global SU(2) gauge transfor- 
mation other than that above changes <5^, so it rotates 
|^o) m to a n e w state orthogonal to \ipo) . The same is true 
for any non-global U(l) gauge transformation (A(f2) not 
constant), because all such transformations change a^. 
We now show that all non-global SU(2) gauge transfor- 
mations change , which implies the result we needed to 
show (i.e. that G either leaves |^o) invariant or rotates 
it into an orthogonal state). 

Under an SU(2) gauge transformation, the object 



transforms as 



(C7) 



In order for a* to 



where U = U(Q) = exp(i\ l (n)a l ) 

be invariant, the above relation must become an equality 
for both [i — 9 and ji = <j) components. Since Ag = 0, we 
therefore must have dglf = 0, and U is independent of 9. 
The /i = 4> equation gives 



1 



2sin( 



t 3 - Ua 3 U^ 



(C8) 



Since the term in brackets on the right-hand is indepen- 
dent of 9, and the left-hand side must be independent of 
9, the only consistent possibility is that both d^U = 0, 
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which implies that Xi is constant and the gauge trans- gauge transformations leaving invariant other than 
formation is global, and also Ua 3 U^ = a 3 , which implies that identified above. 
Ai = A 2 = 0. Therefore, there are indeed no SU(2) 
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